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Abstract. This paper is devoted to the study of quantum dissipation in cluster decay 

phenomena in the frame of the Lindblad approach to quantum open systems. The 

tunneling of a metastable state across a piecewise quadratic potential is envisaged for 

two cases : one and two harmonic wells smoothly joined to an inverted parabola which 

simulates the barrier. The width and depth of the second harmonic oscillator well was 

varied over a wide range of values in order to encompass particular cases of tunneling 

f ■j ■ such as the double well potential and the cluster decay. The evolution of the averages 

^»0 \ and covariances of the quantum sub-system is studied in both under- and overdamped 

^-2 ' regimes. For a gaussian intial wave-packet we compute the tunneling probability for 

different values of the friction coefficient and fixed values of the diffusion coefficients. 
The ansatz used for these coefficients corresponds to the case with temperature T = 
as happens in the cold fission and cluster decay phenomena. 
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1. Introduction 

The question of how collective energy is dissipated in various decaying phenomena 
occuring in solid-state and nuclear physics is a topic of major interest. A particular 
case is represented by the decay of a metastable state through quantum tunneling at 
very small or zero temperatures. Such cold decay phenomena are well known at the 
microscopic level for nuclei undergoing cold fission or emission of heavy clusters |1|, Q] 
but also for macroscopic systems like Josephson junctions || or SQUID's [|J. 

Twenty years ago Caldeira and Legget proposed a model of tunneling of a 
metastable quantum state at T = based on the Feynman path integral technique 
Pj. Assuming a linear dependence on the tunneling coordinate q for the friction 
term entering the Lagrangian, they showed that the dissipation decreases the tunneling 
probability. Since that time a new formalism for studying quatum dissipative systems 
emerged. In this approach the system is opened, i.e. it is regarded as a subsystem 
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interacting with a larger system, referred to as the reservoir or the bath ||. The 
dissipation arises as a consequence of this interaction. 

The Lindblad's axiomatic way of introducing dissipation in quantum mechanics 
consists in replacing the one-parameter group action 

Ut(p) = exp (~Ht) p{t) exp Qift) , (1) 

of a closed physical system with Hamiltonian H, by a new one, $i(p) which describes 
the evolution of the corresponding open subsystem. Like U t , this operator should satisfy 
the conditions of self-adjointness ($^(p) = & t (p)), positivity (3>t(p) > 0) and unitary 
trace (Tr($ t (p)) = 1). Moreover &t(p) —* p when t — > 0. 

Using the above mentioned approach the quantum mechanics of one and two 
coupled harmonic oscillators has been extensively studied || |9|, [HJ. Based on these 
results, very recently, the tunneling probability with dissipation was computed exactly 
for the inverted parabolic potential and approximately for the cubic one |JTT], Q . It was 
found that there might be cases when friction favours the penetration. However the 
inverted parabola is not a good candidate to describe a real tunneling process because 
no metastable state could be conceived in it. The cubic potential is more likely to 
describe the decay of a metastable state but in this case the authors faced the problem 
of dealing with a system of equations in averages and higher momenta which is not 
closed. Therefore all cubic and higher momenta have been droped off providing thus an 
approximate solution. 

In this paper we consider that the quantum decay process, characterized by the 
coordinate q is modelled by the tunneling of a one- dimensional potential barrier whose 
explicite form is given by joining first two parabolic segments as in the Kramers problem 
[13| , and next we add a third parabolic well beyond the barrier. The description of 



tunneling phenomena accompanied by dissipation in Josephson junctions [bl or fission 



of heavy compound nuclei [TIJ offenly employs such barriers. Besides computing the 
averages and covariances of the decay coordinates we are mainly interested to investigate 
the dependence of the tunneling probability on the friction. 

2. Lindblad approach to open one-dimensional quantum systems 

According to the original paper of Lindblad @ the time evolution of an observable 
A(q,p) belonging to a system undergoing a non-equilibrium process is written in the 
Heisenberg picture as follows 

— = l -[H, A\+ l — [q[A,p] + [A,p]q - [A, q]p - p[A, q}) 

- ^[p, [p,A]] - ^[q, [q,A]} + |f ([g, [p, A}} + [p, [q,A]}) (2) 

where A is the friction constant whereas D qq , D pp and D pq are the quantum diffusion 
coefficients. The Hamiltonian H of the one- dimensional subsystem with mass m, 
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corresponding to the above equation of motion, has the general form in coordinate 
q and momentum p 



H=£-+V(q) 

2m 



(3) 



Note that eq.(0) was derived under the assumption of a weak coupling of the subsystem 
with its environment [|16[ . 



Since the dynamics of the subsystem is described in terms of averages and 
covariances, for two self-adjoint operators A and B the following notations are 
introduced : 

a A (t) = Tr(p{t)A), a AB (t) = Tr (p(t) 1 ) - a A {t)a B {t) (4) 

If A = q or p the equations of motion derived from eq. (|2]) looks like 



and 



da q {t) 

dt 
dffp(t) 

dt 

do- qq {t) 
dt 

d(Tpp(t) 
dt 



d(Tp q (t) 
dt 



XaJt) H a p (t) 



m 



- Tr (" w ^)-^ (i) 



- 2\a qq (t) + -a pq (t) + 2D qq 

- 2Xa pp (t) - Tr L(t)^^-p + h.c\ 

+ 2Tr(p(t)^^\a p (t) + 2D pp 

~ 2Xa„(t) + ±cr„{t) + a q (t) Tr (p(t)^^j 

-TrUt) d -^±q) + 2D pq 



(5) 



(6) 



The above set of first-order o.d.e. has been solved exactly for the harmonic oscillator, 
V = mu 2 q 2 /2, and the inverted parabolic potential (u — > io;)|| [12 . 



3. Quantum tunneling across quadratic potentials 



We formulate the decay problem as the tunneling of an initial gaussian wave packet, 
confined in a harmonic oscillator well (potential pocket), located at the left of the barrier, 
with given values of the position average c q (0) and covariance c qq (0) : 



ijj{q) 



(2tt^(0))V4 



exp 



4M0) 



(q - a q (0)) 2 + -a p (0)q 



n 



(r) 



According to |17|] , the transition probability for an irreversible process described by 
the set (||,|D reads 



W(q,p,t) = (2vr)- 1 (Det a 



v-l/2 



exp<j-- 



\t)]jq-o- q (t)f 
1 



°- -1 (*)j M (? - ^(t))(p - op®) - 2 r^\ PP (p - ^yy 



(8) 
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where cr is the 2x2 matrix of the covariances defined in eq. (HI) . The probability to find 



the wave-packet to the right of the barrier is [18 



P{q b -t) = f +0 °dq f +C °dp W(q,p,t) = ^erfc fa® =B ) (9) 



The decay rate, or fission width, is given by the ratio between the diffusion current 
across the barrier, J(q b ; t) = f_™ dp W(q b ,p; t) and the tunneling probability 

J {(lb] t) 



r/(t) 



P(g b ;t) 

1 ex p 



KW-gh) 2 



2a qq (t) 



, a qq (t)a p (t) + a pq (t)(q b - a q (t))) L ^ J (10) 

2na qq (ty erfc ( ^kSL ) 

The purpose of the next two subsections is to develope the above formalism for two 
case studies, namely the tunneling of a metastable state from a parabolic pocket to a 
quadratic unbound potential and next the tunneling between two quadratic wells with 
finite depths, separated by an inverted parabola. 

3.1. Two smoothly joined parabolic potentials 

In this case the potential is given by joining at q t , an inverted oscillator to an upright 
ground-state oscillator (see Fig.l): 

V{q) = V(q a )+ l -n\{q-q a )\ q < q t 

= V(q b )-±n 2 b (q b -q) 2 , q > q t (11) 

In the above expression one suppose that the height of the barrier B = V(q b ) — V(q a ), 
and the frequency Q b of the inverted oscillator are known. In the case of nuclear reactions 
this means to know the barrier from the heavy-ion interaction of the two clusters or the 
dependence of the total energy on deformation |19|, [2(J. Assigning a decay energy E to 



the particle of mass m we suppose that the minimum of the first harmonic oscillator well 
q a coincides with the first turning point. Naturally, since we deal with smoothly joined 
segments, the continuity of V and its derivative dV/dq is required at q = q t . From here 
we get the frequency of the well's bottom : 

a = nl(q b -q a y-2B (12) 

and the position of the joining point : 

q a Ql + q b n% 
qt - w a + ni (13) 



In order to integrate the systems (|^J§) we need to establish the initial conditions. 
The initial wave-packet will be confined inside the harmonic oscillator well centered 
at o"g(0). In our calculation, motivated by specific examples from nuclear fission we 
choosed cr g (0) = q a) i.e. the gaussian wave-packet is centered at the minimum of the 



Decay with dissipation 5 

potential pocket. In order to overcome the barrier the initial momentum should satisfy 
the inequality cr p (0) > \JlmB for a vanishing friction coefficient A. The initial value 
of the coordinate covariance o~ qq (0) is choosen in such a way that its time derivative 
vanishes. Assuming that o~ pq (0) vanishes too, we then get from eq.(^) 

MO) = ^ (14) 

Since we deal with a gaussian wave packet, the momentum and the coordinate 



covariances are simply related by means of the following relation [T2j 



ff " (0) = sSii) (15) 

For the diffusion coefficients we take the rotating-wave approximation [ EH] 

Xh 

D qq = — -^=D pp = mQ 2 a D qq , D pq = (16) 

2 A raVtl 



These values were obtained from eq.(|B|) by asuming that the initial Gaussian metastable 
state is a Gibbs state at T = 0. 

In Figs. 2 and 3 we represented the averages and covariances of the coordinate and 
momentum taking four choices of the friction coefficient A. We assumed that the initial 
momentum is the same, i.e. <t p (0) = 1200 MeV. The general trend is a damping of the 
motion due to friction. Above a critical value A cr of the friction coefficient all the 
momenta tends to infinity, faster if A = and slower if A /* A cr . Bellow this critical 
value the gaussian wave packet will perform damped oscillations around the potential 
minimum. From the inspection of Figure 2 it is obvious that the wave-packet is spending 
a longer time in the barrier if A is approaching its critical value. The decreasing of the 
tunneling probability (see eq.(j^)) with A is shewn in Fig. 4. After a certain transition 
time (« 5 x 10 _22 s) the tunneling probability will reach a non- vanishing asymptotic value 
for A < A cr . Looking on the first panel of Fig. 2 one sees that after this lapse of time 
the centroid of the gaussian is located on the other side of the barrier. Consequently, in 
the present approach, the tunneling time, a different quantity from the metastable state 
life-time, is assigned to the time necessary to attain an asymptotic value P(qb] oo). For 
a discussion of this subject in the frame of the WKB-approximation and time dependent 
approach to Schrodinger equation, see [f22], f2"3"f . 

For A > A cr the tunneling probability eventually tends to zero but at the beginning 
of the motion, between certain lapses of time, the right tail of the wave function is 
found inside the barrier. If the initial momentum cr p (0) would be increased then the 
barrier is once again overcomed. However in tunneling phenomena like cluster decay 
the initial average momentun of the wavefunction should not exceed certain limits and 
thus upon comparison with the penetrabilities deduced from experiment, the range of 
possible values of the friction coefficient could be deduced. 

Thus, for two smoothly joined parabolas, the friction reduces the tunnelling 
probability of the metastable state in the frame of Lindblad's axiomatic approach. 
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3.2. Three smoothly joined parabolic potential 

The case studied bellow suffers from the deficiency that the right tail of the potential 
falls to — oo. In usual applications we deal with tails which tends to +oo or to 0. We 
thus join a third parabolic segment at a second turning point, q t2 , centered at q c (see 
fig.5) : 

V{q) = V{q a ) + hl 2 a {q-q a )\ q < q tl 

= V{q h )-hll(q h -q)\ q tl < q < q t2 (17) 

= V(q c ) + -n 2 c (q-q c ) 2 , q > q t2 

Imposing the same conditon of continuity of V and its derivative dV/dq at q = qt 2 one 
gets 



ni 



& b (q c -q b y-2AV l 



be 



and 



_ g c n 2 c + qb n 2 b 
qt2 - ni + ni (19) 

where /S.V bc = V(qj,) — V(q c ) is the difference between the top of the barrier and the 
bottom of the second well. 

In Fig. 6 the evolution of averages were represented for two harmonic wells with 
equal bottoms, i.e. V(qb) = V(q c ) (a) and V(q c ) = (b). For the location of the second 
well minima several values were selected. The first one (q c = 16.5 fm) corresponds 
almost to the symmetric double-well problem. In the absebce of friction the tendency 
of o~ q (t) is to oscillate indefinetely between these two wells with amplitude and period 
determined by the location of the second minimum q c . Like in the case examined in the 
preceeding subsection, a critical value of A will occur. Above this value the wave-packet 
will never cross the barrier and will perform damped vibrations inside the first well. If 
A < A cr then the gaussian will pass in the second well and will not return in the first 
well due to the kinetic energy loose by friction. In this last case the centroid of the 
wave packet will tend asymptotically to the center of the second well q c . Ase one sees 
from Fig. 6 there are no qualitative difference between the motions in the potential with 
V(q b ) = V(q c ) and V(q c ) = 0. 

In terms of tunneling probability, the results obtained for the three smoothly joined 
parabolas are displayed in Fig. 7. If friction is absent the tunneling probability will tend 
to 0.5, i.e. half of the wave-packet is found at the right of the barrier's top, q b , and half 
to its left. If friction is switched on, but up to A cr , then asymptotically the wave-packet 
tunnels 100% in the second well after a certain lapse of time. This lapse of time will be 
shortened if the minima of the second well is shifted towards higher values. When A cr 
is exceeded the tunneling probability goes to zero since the wave packet is constrained 
by the high friction to not leave the first harmonic well. 
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Contrary to the preceeding case, for the double- well problem the value of the friction 
constant, bellow a threshold, favours the tunneling. Above this threshold the tunneling 
is totally hindered. 

4. Conclusions 

In this paper the study of tunneling, in a one-dimensional piecewise potential built-up 
from two and three parabolic segments, was carried out in the frame of the Lindblad 
approach to quantum systems for a gaussian wave-packet. In this way the interaction 
of the system with the environment was taken into account and the effect of dissipation 
incorporated in the equation of motion for the averages and the covariances. In view of 
the applications in cold fission and cluster decay only the case with temperature T=0 
was considered. This amounts to admit a linear dependence of the diffusion coefficients 
on the friction constant for an initial metastable state. 

For the two-smoothly joined parabolic potentials, if the friction constant is bellow 
a certain critical value, the gaussian wave-packet cross just one time the barrier and 
afterwards it moves to infinity with a high or a low momentum, depending on how large 
is the friction. The computation of the tunneling probability in this case confirms older 
results || which claimed that the effect of linear friction decreases the chance to find 
the wave-packet beyond the barrier. 

In the case of three-smoothly joined parabolas the situation is different. If the 
motion is unviscous then the wave packet will oscillate unhindered between the infinite 
walls of the two wells. Therefore, the probability current flowing from the left to the 
right of the barrier will be compensated by the one flowing in the opposite direction. 
If the friction is present in the system, although bellow a certain limit, the tunneling 
is stimulated, in the sense that asymptotically all the wave-function is located in the 
right well. This happens no matter how large or deep is this second well. Only the 
time necessary to establish the equilibrium will depend on the features of the employed 
potential. 

It would be interesting to extend the calculations presented in this paper to other 
one-dimensional potentials which are more likely to simulate the behaviour of systems 
undergoing decay by quantum tunneling. One of the limitations of the present approach 
is that the potential needs to be provided in analytical form, preferably as a polynomial 
of lower degree. Already for a cubic polynomial complications occur because the system 
of first order differential equations describing the time evolution of averages and higher 
momenta is non-linear and is not closed. A solution would be to consider the realistic 
one-dimensional potential as a smooth piecewise curve given by cubic splines. Otherwise, 
the exact treatment of the problem would be to solve directly the time-dependent 
Schrodinger equation, which unfortunately is non-local and non-linear [Rf. 
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Figure 1. Gaussian wave-packet with average initial momentum <r p (0) centered at q a 
undergoing tunneling across a barrier of height B built-up from two smoothly joined 
parabolas. 



Figure 2. The time dependence of the average values of the coordinate a q (t) and 
momentum cr p (t) for two smoothly joined parabolas. Four different values of the friction 
constant A are taken. 



Figure 3. The evolution of <r qq (t), cr pp (t) and cr pq (t) for friction constants taken as in 
Fig.2. 



Figure 4. The decrease of the tunneling probability in time with the friction constant 
in the two smoothly joined parabolas case. 



Figure 5. Same as in Fig.l but for three smoothly joined parabolas. The location q c 
and depth of the second well minima arc varied over a wide range. 



Figure 6. The evolution of cr q (t) for three smoothly joined parabolas. Three different 
values are taken for the friction constant A and four choices for q c . In (a) the depth of 
the two harmonic wells are equal, i.e. V(q a ) = V(<?&), whereas in (b) V{q c ) = 0. 



Figure 7. The tunneling probability in the three smoothly joined parabolas case for 
(a) V(q a ) = V(q b ) and (b) V(q c ) = 0. 
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